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ABSTRACT 

Observations of early- type galaxies, both in the local universe and in clusters at 
CN ■ medium redshifts, suggest that these galaxies often contain discs or disc-like struc- 

tures. Using the results of Kelson et al. (2000) for the incidence of disc-components 
among the galaxies in the redshift z = 0.33 cluster CL 1358+62, we investigate the 
effect of disc structures on the lensing properties of early- type galaxies. Statistical 
properties, like magnification cross sections and the expected number of quad (four- 
image) lens systems, are not affected greatly by the inclusion of discs that contain less 
| than ~ 10 per cent of the total stellar mass. However, the properties of individual lens 

f^i . systems are affected. We estimate that 10 — 30 per cent of all quad lens systems, with 

early-type deflector galaxies, would be affected measurably by the presence of disc 
components. Intriguingly, the image magnification ratios are altered significantly. The 
amplitude of the predicted change is sufficient to explain the observed magnification 
ratios in systems like B1422+231 without requiring compact substructure. Further- 
more, time delays between images also change; fitting a bulge-only model to early-type 
lenses that in fact contain a disc would yield a value of the Hubble constant Hq that 
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is systematically too low by ~ 25 per cent. 
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1 INTRODUCTION 

In current models of hierarchical galaxy formation the mor- 
phological type of a galaxy is determined to a large extent by 
its merger history (Abraham & van den Bergh 2001): late- 
type galaxies have not undergone a recent massive merger, 
so that stable, prominent discs formed, whereas early-type 
galaxies have been prevented from forming a prominent disc 
by a major merger that either consumes all the available 
gas for star formation, or redistributes angular momentum 
via large scale winds (Somerville & Primack 1999). This be- 
haviour has also been recovered in hydrodynamical simu- 
lations of galaxy formation (Steinmetz & Miiller 1995). So 
far, there is no evidence for any evolution in the luminos- 
ity and size of discs in late-type galaxies up to redshifts 
of about 1 (Simard et al. 1999; Ferguson & Clarke 2001). 
This suggests that if mergers occur after that time, they 
are either not massive enough to significantly affect the disc 
structures or they destroy visible disc structures, producing 
elliptical systems. It is not known on what time scales disc 
structures form. Simulations and some observations suggest 
that disc structures form slowly between redshifts 0.1 and 



1 (Mo, Mao & White 1998) but there is also evidence that 
some large discs are already present at higher redshifts (van 
Dokkum & Stanford 2001). Most likely, discs in early-type 
galaxies are the result of remaining gas settling after ma- 
jor merging events between two late-type galaxies (Naab & 
Burkert 2001). Another possible formation scenario for discs 
involves slow accretion from many minor merging events. 
Debris from minor mergers, or high velocity clouds in the 
halo of the galaxy, could rain down onto the disc (Benjamin 
& Danly 1997; Putman et al. 2002) . However, such a scenario 
would lead to thicker, less pronounced disc structures. Since 
late-type galaxies exhibit clear disc structures, this process 
is unlikely to be very important in the formation of typical 
spiral structures. 

In all these scenarios, the obvious question arises as to 
whether disc structures are exclusively found in late-type 
galaxies. Observationally, that is clearly not the case. SO 
galaxies in the local universe possess disc components, and 
even local ellipticals show evidence of disc structures and 
dust lanes when studied in detail (Rest et al. 2001; Tran 
et al. 2001). The isophotes of ellipticals are usually classified 
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into "boxy" and "discy", types. Pointed, discy isophotes may 
well be caused by the presence of discs. Even in the case 
when the isophotes are nearly round and show no sign of 
"discyness" , kinematic studies suggest that it is still possible 
to hide a disc in such galaxies (Romanowsky & Kochanek 
1997). 

There is also some observational evidence for discs in el- 
lipticals at cosmological redshifts; Kelson et al. (2000) found 
that many early-type galaxies in CL 1358+62 at z = 0.33 
possess evidence for the presence of exponential discs, sug- 
gesting discs are present in a large fraction of early-type 
galaxies, at least in clusters. Since there is no clear under- 
standing yet on how these discs form, and no systematic ob- 
servational investigation into the abundance and properties 
of discs in early- type field galaxies has been carried out, it is 
difficult to make precise statements about the abundance of 
discs in field ellipticals. However, some recent observations 
using the Hubble Space Telescope (HST) show that a frac- 
tion of at least ~ 50 per cent of field ellipticals contain discs 
of sizes spanning a range of at least an order of magnitude, 
from small nuclear discs to large discs extending beyond the 
bulge (Rest et al. 2001). In addition, kinematic studies using 
the SAURON integral field spectrograph, also suggest that 
a high fraction of ~ 50 per cent of field ellipticals contain 
discs of various sizes (de Zeeuw et al. 2002) . 

Most lens galaxies known to date are early-type galax- 
ies. The main reason for this is that early-type galaxies tend 
to be the most massive and hence have a larger cross sec- 
tion for gravitational lensing. However, as shown by Blain, 
Moller & Mailer (1999) spiral galaxy lenses are also pre- 
dicted to contribute significantly to the total lensing cross 
section. This is mainly due to the high surface mass density 
of the thin disc component when viewed edge on (Mailer, 
Flores & Primack, 1997; Moller & Blain, 1998). As shown 
in those papers, the disc components of late-type lenses can 
make important contributions to the lensing potential and 
change the generic properties of the lens systems. The dis- 
covery of spiral lens systems in the CLASS lensing survey, 
e.g. B1600+434 (Koopmans, de Bruyn & Jackson 1998), 
confirm that this is the case. However, early-type lensing 
galaxies have so far nearly always been modelled with a sim- 
ple one-component model, possibly including external shear 
(Keeton, Kochanek & Seljak, 1997; Lehar et al., 2000). 

In this paper we determine the lensing properties of 
early-type galaxies that contain realistic discs with stellar 
masses of up to ~ 10 per cent of that of the bulge. We 
will determine the lensing properties of early-type lenses and 
show how discs in ellipticals could explain some properties 
of observed lens systems that have not yet been modelled 
successfully. We concentrate on three main questions: 

(i) statistical lensing properties: how do disc structures 
in early-type lensing galaxies affect the statistical proper- 
ties of a sample of early-type lenses? We calculate magni- 
fication cross sections and the expected ratio of quad: two- 
image systems. In particular, we investigate whether discs in 
early-type lenses may explain the high ratio of quad image 
systems, which is still an unsolved enigma in gravitational 
lensing (Rusin & Tegmark, 2001, but see Chae, 2002). 

(ii) magnification ratios: can disc structures explain some 
or all of the unusual magnification ratios observed in some 
early-type lens systems? To date, these deviant magnifica- 



tion ratios have been interpreted as providing strong evi- 
dence in favour of the presence of compact halo substructure 
as predicted by current cold dark matter (CDM) simulations 
(e.g. Metcalf & Madau 2001). 

(iii) time delays: how would disc structures affect the 
measured time delays in early-type lens systems? Currently 
there seems to be a strong disagreement between the value 
of Ho as found in the HST Key Project (Freedman et al. 
2001) and measurements of Ho from gravitational lensing 
(Kochanek 2003). Some local estimates of Ho (e.g.Saha et al 
2001) also favour lower values. We investigate whether disc 
structures in gravitational lens systems may reduce this dis- 
crepancy. 

We start in §2 by discussing the appearance of the model lens 
galaxies and show that the manifestations of the presence 
of a disc are similar to those observed in some ellipticals. 
In §3 we present the parametric lens models used in this 
paper. In §4 we describe some general and statistical lensing 
properties of early-type galaxies that include disc structures. 
In §5 we show how discs in ellipticals could explain some of 
the observed magnification ratios and in §6 we calculate how 
the presence of discs in ellipticals would affect estimates of 
Ho from lensing. We conclude with a summary in §7. 

Throughout this paper we assume a Friedmann- 
Robertson- Walker Universe with Om = 0.3, £!a = 0.7 and 
Ho = 50kms~ 1 Mpc -1 . The choice for this value of Ho is 
motivated due to limitations of the original lensing code. 
Subsequently, the code has been extended to allow for any 
value of Ho, as necessary for the calculations in §6, but we 
did not recalculate the results of §3-5, as these calculations 
are very time-consuming and a change in Ho has little effect 
on the results. A value of Ho = 70kms _1 Mpc -1 , is prob- 
ably more in line with current observational constraints, as 
given for example by the HST Key Project (Freedman et al. 
2001) and the WMAP results (Bennett et al. 2003). 



2 SURFACE BRIGHTNESS PROFILES 

The existence of disc components with large scale-lengths 
in luminous elliptical galaxies has been known for some 
time (Bender et al. 1989; Rix & White 1990; Rix 1991). 
Recent theoretical work (Naab & Burkert 2001; Barnes 
2002) has begun to address the origin of such components 
in bulge-dominated galaxies. However, the ability to reli- 
ably determine the prevalence of discs in strongly bulge- 
dominated systems has only been possible following the de- 
velopment of full two-dimensional decomposition techniques 
(Byun & Freeman 1995). Given data of adequate signal-to- 
noise ratio and angular resolution, two-dimensional decom- 
position methods are capable of identifying discs in systems 
with bulge fractions (B/T) as high as B/T=0.95. A signif- 
icant quantity of HST imaging, obtained using the Wide 
Field Planetary Camera 2 (WFPC2), of luminous bulge- 
dominated galaxies in intermediate redshift clusters is well- 
suited to such decompositions but little work appears to 
have been undertaken. Exceptions are the investigations of 
the z = 0.33 cluster CL1358+62 by Kelson et al. (2000) and 
Tran et al. (2003). Fig. 4 of Kelson et al. provides a partic- 
ularly striking demonstration of the presence of discs, even 
among galaxies with conventional 'E' morphological classi- 
fications. 
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In order to parameterize the light profile of ellipticals, 
a commonly used set of parameters are the half-light radius 
rh and the surface brightness at the half light radius Ih- 
The surface brightness of a de-Vaucouleurs bulge can then 
be written as (Binney & Tremaine 1987), 

I(r) = I h cxp{-7.67[(r/r h ) (1/4) -1]}. (1) 

The parameterization of the surface brightness distribution 
for an exponential disc, in terms of the central surface bright- 
ness, 7d, and scale-length, ra, can be written 

I{r) =/ d exp-(r/r d ). (2) 

The bulge fraction, B/T, for a galaxy parameterized by 
such bulge and disc components is then 



T I h rl+0.28I d r 2 d w 

where Ih and rh, Id and r<j are as defined above. 

Table 2 in Kelson et al. gives surface brightnesses and 
half-light-radii for bulge-only and bulge-plus-disc fits to 
galaxies in CL1358+62. Note that Kelson et al. specify the 
properties of the disc component in terms of half-light radius 
of the disc, rhd = 1.688rd, and the mean surface brightness 
within the half-light radius, < Jhd >• Using this parameter- 
ization the B/T ratio may be calculated by replacing the 
term 0.28/drd by < J hd > r'l d . 

For the investigation of the general and statistical lens- 
ing properties of discy ellipticals, we use the full sample 
of galaxy models in Kelson's catalogue. However, our main 
results, as presented in § 5 and § 6 involve time-consuming 
calculations and are therefore based on the parameters de- 
scribing a single model galaxy. The presence of many galax- 
ies with bulge and disc components with comparable scale- 
lengths in CL1358+62 is clear from the form of the residuals 
shown in Figure 4 of Kelson et al. It is worth stressing that 
the presence of even highly inclined (almost edge-on) discs, 
with scale- lengths similar to those identified by Kelson et al., 
making up less than ~ 10 per cent of the total light are hard 
to discern directly in images. Plots showing the residuals ob- 
tained by fitting a bulge-only model are, however, far more 
revealing. The presence of face-on discs is revealed by char- 
acteristic donut-like residuals, while edge-on discs appear as 
often quite prominent inclined structures. Notwithstanding 
the evidence provided by Kelson et al. in support of the ac- 
curacy of their bulge-disc compositions we have chosen to 
be conservative by reducing the disc contribution to each 
galaxy in the Kelson et al. catalogue by half. Thus, our in- 
vestigation of the statistical lensing properties is based on 
the properties of the galaxies listed in Table 2 of Kelson et 
al. but with a reduction in the disc contribution of a factor 2. 
For the individual model galaxy used in our calculations we 
have chosen galaxy no. 242 (Table 1). Adopting our rescal- 
ing of the disc contribution, galaxy no. 242 has B/T=0.95, 
c.f, B/T=0.91 in Kelson et al. Galaxy no. 242 has a total 
mass, B/T ratio and bulge scale-length close to the average 
value for the 'E' galaxies in Kelson et al's catalogue. In our 
simulations the ellipticity of the bulge is set to e = 0.4 and 
the inclination of the disc is i — 75 deg. 

Fig. 1 shows images based on our adopted parameteri- 
zation of galaxy no. 242 with B/T=0.95. Fig. la is an image 
based on the bulge-only model surface-brightness profile fit 
(parameters from Cols: 3 & 4 in Table 1), Fig. lb shows a 



realization of the bulge-plus-disc surface-brightness profile 
(parameters from Cols: 7 & 8 (bulge) and Cols: 10 & 11 
(disc) in Table 1). The disc is at high inclination, i = 70 deg 
to the line of sight (i — 90 deg corresponding to edge-on) 
and oriented at an angle of 45 deg with respect to the x- 
axis. The presence of the disc is barely discernible. Fig. lc 
shows the residual image derived from subtracting Fig. la 
from Fig. lb (both normalized to possess the same total 
counts). The images may be compared to those in Figs 1 
and 4 of Kelson et al.. In Kelson et al., galaxy no. 242 shows 
the donut-like residual characteristic of the presence of a 
close to face-on disc. By contrast, in Fig. lc, the disc ap- 
pears prominently because of the high inclination angle of 
i — 70 deg to the line-of-sight. 

Following submission of the original version of this pa- 
per an analysis of the photometric properties of galaxies in 
CL 1358+62 undertaken by Tran et al. (2003) appeared. 
Comparison of both the B/T values and the scale-lengths 
of the light distributions for the bulge-dominated galaxies 
in common between Kelson et al. and Tran et al. show a 
good correlation over the full range of B/T. However, for 
large values of B/T the contribution of the disc components 
are much larger in the Tran et al. analysis. Less satisfac- 
tory still is the comparison between the scale-lengths of the 
disc components found among the bulge-dominated galaxies. 
The disc scale-lengths found by Tran et al. are significantly 
larger, explaining, at least in part, the systematically smaller 
B/T ratios found for the bulge-dominated systems. Tran et 
al. present the results of a large number of Monte-Carlo sim- 
ulations designed to quantify the accuracy of the decompo- 
sitions into bulge and disc components. No reference is made 
to the strong systematic disagreement between the statistics 
of the disc scale-lengths found in the two studies. Tran et al. 
are confident that their simulations demonstrate the ability 
of the software employed to obtain accurate bulge-disc de- 
compositions. However, the parameter ranges employed in 
the simulations appear to be determined from the results of 
the actual decompositions applied to the real data. Tran et 
al. find no discs with scale-lengths comparable to or smaller 
than those of the bulge components of the galaxies. Thus, it 
is not clear that the ability of Tran et al.'s analysis to cor- 
rectly parameterize combinations of bulge and disc profiles 
with parameters found by Kelson et al. has been adequately 
tested. In the circumstances we have decided to retain the 
very conservative distribution of B/T ratios, with associated 
values of component scale-lengths, based on the distribution 
advocated by Kelson et al. 



3 THE LENS MODELS 

In order to study the lensing properties of early-type galax- 
ies including discs, a working model of the mass distribution 
is needed. Unfortunately, little is known about the detailed 
properties of early-type galaxy discs and how they depend on 
morphology and environment. From past studies it is clear 
that discs in ellipticals are abundant and that the disc sizes 
and magnitudes cover the whole range from compact and 
luminous to extended and faint (Rest et al. 2001; Scorza & 
Bender 1995; Scorza et al. 1998). Most of the available infor- 
mation is based on observations of individual galaxies; there 
has not yet been a systematic high-resolution study of the 



4 Ole Moller, P. Hewett & A.W. Blain 



7 



1 1 1 1 1 F 1 1 1 1 1 1 1 1 1 1 1 1 1 

: (a) : 


.iii 

: (b> 


— i — i — i — i — | — i — i — i — i — | — r — i — i — i — 



















_"L 


Bulge only - 
i i i i l i i i i l i i i i l i i i i 


i i i i 


Bulge+disc j 




x / arcsec 



x / arcsec 



x / arcsec 



Figure 1. Simulated images of a typical model galaxy at redshift z = 0.33. The galaxy model parameters correspond to our modified 
representation of galaxy no. 242 in Kelson's catalogue (Table 1). Panel (a) shows the image derived using a bulge-only model. Panel (b) 
shows the bulgc-plus-disc model, normalized to have the same total counts as the bulge-only model. The disc component of the model is 
inclined at 70 deg to the line of sight and oriented at an angle of 45 deg, to the x axis. Panel (c) shows the residual following subtraction 
of the bulge-only image from the bulgc-plus-disc image. The pixel scale, 0.1 arcsec per pixel, and point spread function have been chosen 
to approximate images obtained using the WFPC2 on HST. The signal-to-noise ratio of the synthetic images is significantly higher than 
that typical of HST imaging data of the type analysed by Kelson et al. (2000). 



morphologies of a complete, unbiased sample of early-type 
galaxies. In order to make meaningful predictions about the 
lensing properties of early-type lensing galaxies including 
disc components, we base our models on the light profiles 
obtained in Kelson et al. (2000). All galaxies in Kelson's 
sample are at a redshift of z = 0.33, which corresponds 
roughly to the redshift at which lensing is most efficient; 
few lenses are expected with redshifts less than z < 0.1. 

3.1 Parameterization of early-type lenses with 
discs 

We model the elliptical lensing galaxies with two- 
components: an elliptical bulge with a de-Vaucouleurs r 1 ^ 4 
profile and an exponential thin disc. 

Under the usual assumption of a thin lens, the lensing 
properties depend solely on the surface mass density profile. 
For the bulge component this is given by 

E b (r) =S b exp(-r7r b ) 1/4 , (4) 

where E b is the central surface mass density of the bulge, 
r b is the bulge scale-length, and 

r = y/x 2 + y 2 (l -e) 2 . (5) 

Here r = x 2 + y 2 , and c is the ellipticity of the bulge. 
The disc, inclined at an angle i to the line of sight, can be 
parameterized in a similar way, to give 

S d (r) = S d exp (-r'/r d ), (6) 

where Ed is the central surface mass density of the disc, rd 
is the scale-length of the disc and 

r = ^/:r 2 + j, 2 [l-cos (i)} 2 . (7) 

We convert the light profiles as given in equations (1) & 
(2) to models of the mass distribution assuming a constant 
mass-to-light ratio, T. The relation between the parameters 
in equations (1) & (2) and those in equations (4) & (6), is 



then r h = 7.676 4 r b and J h = r _1 E b exp (-7.676) for the 
bulge component and Id = r _1 £d for the disc component. 
We will refer to rh and rhd = 1.688r d as the 'effective ra- 
dius' or 'half-light radius' of bulge and disc, respectively. 
Assuming randomly oriented discs, the average inclination 
is about 70 deg. Unless otherwise stated, the discs in our 
bulge-plus-disc models have this inclination. 



3.2 Possible caveats of the model 

The photometric profile decompositions of galaxies in CL 
1358+62 provide an observational motivation for mod- 
elling early-type galaxies including disc structures. However, 
galaxies in cluster environments might exhibit more pro- 
nounced discs than corresponding galaxies in the field. If 
this is indeed the case, and simulations of Naab & Burk- 
ert (2001) suggest otherwise, the discs in our model lensing 
galaxies could be too pronounced compared to those ex- 
pected in observed early-type lens systems, many of which 
are in lower-density field environments. A comparison of the 
disc sizes obtained from a sample of mostly field ellipticals 
by Scorza & Bender (1995) with the sizes of discs in E's 
and SO's found in Kelson et al. (2000) shows that the discs 
from both samples are similar and range in disc half-light 
radii from 0.1 kpc to lkpc. This could indicate that the en- 
vironment may influence the relative number of galaxies as 
a function of morphology, but has little effect on the prop- 
erties of discs in galaxies of a given morphological type. In 
any case, we are not concerned here with accurate modelling 
of the statistical incidence of discs, but rather aim to show 
how discs with sizes and masses consistent with current ob- 
servations influence the gravitational lensing properties of 
early- type lens systems. 

There is a second possible caveat with the mass model 
we use here: we assume that mass follows light and do not 
include a dark matter halo. As we will discuss in more de- 
tail below, the main effect of inclined, thin disc components 
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Figure 2. Radial surface mass density profiles for representative 
models of the galaxies in the cluster CL 1358+62. The different 
lines show the surface mass density as a function of radius for dif- 
ferent model parameters that are representative for the majority 
of all galaxies in the sample. All curves are normalized so that 
the total mass is lO n M . The ellipticity of the bulge is e = 0.4 
and the inclination angle of the disc is i = 70 deg. The surface 
mass density differs by less than 10 per cent between models that 
include a disc with M<j = 0.1Mb and models without a disc. At 
small radii, discs increase the surface mass density, due to their 
smaller scale lengths. At the cluster redshift of z = 0.33, 1 arcsec 
corresponds to 6.6 kpc for the cosmology we adopt in this paper. 

on the lensing properties of galaxies arises from the strong 
asymmetry that is introduced in the central regions of the 
galaxies due to the disc. We therefore expect that the pre- 
cise form of the spherical or nearly-spherical components in 
which the disc is embedded, that is, the bulge and halo, will 
not change the effect of the added disc component drasti- 
cally. We have checked this explicitly for the mass model 
used to derive the results in §5 and §6, by adding a dark 
matter halo that contains 2/3 of the total mass and scaling 
the mass of the bulge and disc so that the Einstein radius 
is the same as for the model without a dark-matter halo. 
Estimates of the mass-to-light ratio of early type galaxies 
range from M/L ~ 2 in the K-band (e.g. Moriondo, Gio- 
vanardi & Hunt 1998) to about M/L ~ 10 in the V-band 
(Natarajan et al. 1998). Many lens systems are consistent 
with lower mass-to-light ratios within an Einstein radius of 
about M/L ~ 2 in V (Koopmans & Treu 2003), motivating 
our choice of a dark halo containing 2/3 of the total mass. 
When including such a spherical dark halo we find that the 
inner caustic shrinks in size, and so quad images become less 
likely, but that the magnification ratios and time delays of 
images of sources that are quadruply imaged do not change 
substantially. 

3.3 Surface density and mass profiles of lens 
models 

We show the radial surface mass density profiles for some 
typical cases in Fig. 2. The solid lines show the density pro- 
files of three bulge-only models with effective radii of 0.2, 1.0 
and 2.0 arcsec. At the cluster redshift of z = 0.33 these val- 
ues correspond to 1.3, 6.6 and 13.3 kpc respectively. These 
values cover the range of half-light radii of the bulge-only fits 
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Figure 3. Radial mass profiles for a number of galaxy models. 
The different lines show the enclosed mass as a function of radius 
for different model parameters. The curves are normalized so that 
liniR—xx; M(R) = 1O 11 M0. The effective radii for bulge and disc 
are rhd = Hi = 0.2 arcsec, which corresponds to 1.3 kpc at the 
cluster redshift of z = 0.33. The ellipticity of the bulge is e = 
0.4 and the inclination angle of the disc is 75 deg. The enclosed 
mass at distances of ~ 1 arcsec differs by less than 10% between 
models that include a disc containing 10% of the mass and models 
without a disc. Note that, as in the previous figure, the most 
significant contribution of the disc is in the inner regions. 

to the light profiles of the galaxies in CL 1358+62, which 
lie between 0.2 and 4.9 arcsec with only two galaxies hav- 
ing half-light radii larger than 1.6 arcsec. The dashed lines 
in Fig. 2 show the surface mass densities for discs of effec- 
tive radii of 0.05, 0.2 and 0.3 arcsec. The discs are assumed 
to be inclined at 70 deg to the line of sight. Most of the 
bulge-plus-disc fits to the light profile of the elliptical galax- 
ies in the sample have discs with effective radii less than 
0.3 arcsec, corresponding to 2 kpc at z = 0.33. The disc com- 
ponents of the later type galaxies are usually larger, with 
effective radii of up to 1 arcsec. Comparison of the relative 
bulge and disc contributions shows that the disc components 
contribute significantly to the surface mass density within 
about 0.5 arcsec, or 3.3 kpc. The dotted lines show the to- 
tal surface mass density for a typical bulge-plus-disc fit, but 
with different B/T mass ratios. It is clear that the radial 
profiles for models that include a disc of less than ~ 5 per 
cent of the total mass are very similar to bulge-only mod- 
els. For B/T mass ratios of 0.9 or less, that is for most of 
the SO/Sa or later type galaxies in Kelson et al., the surface 
mass density in the inner 0.2 arcsec differs by more than 
10 per cent. 

Gravitational lensing measures directly the projected 
mass Menc contained inside a circular region of radius R. 
In Fig. 3 we show M cnc as a function of radius R. All the 
models shown are normalized to the same total mass, but 
the B/T mass ratio varies. The effective radius of disc and 
bulge, ellipticity of bulge and inclination angle of disc are 
Hid = 0.2 arcsec, r h = 0.2 arcsec, e = 0.4 and i = 75 deg 
respectively and are the same for all models. These values 
correspond to the typical values for an elliptical galaxy in CL 
1358+62. At distances of R ~ 1 arcsec the enclosed mass for 
typical bulge-plus-disc models with B/T> 0.9 is only a few 
percent larger than the enclosed mass for bulge-only models. 
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None of the discs in the light profile fits to the early-type 
galaxies in CL f358+62 contain more than ~ 10 per cent of 
the light (cf. Table Al), and so, if the mass-to-light ratio 
is constant, then disc structures are expected to contain at 
most 10 per cent of the mass. The effect of the disc on the 
maximum image separations or the radius of the Einstein 
ring will be small in those cases. More massive discs with a 
mass of Md ~ 0.3Mb, contribute significantly to the mass 
in the inner region of the galaxy and influence the mass 
contained within the expected Einstein radius more strongly. 

The properties of all the elliptical galaxies from the Kel- 
son et al. (2000) catalogue with the reduced disc contribu- 
tion adopted for our simulations are summarized in Table 1. 
The table lists both the values for the best fitting bulge-only 
models and the values for the best bulge-plus-disc fits. The 
two parameters f c and ft give the B/T mass ratio as mea- 
sured at the Einstein radius _Re and at R — > oo respectively. 
The Einstein radius is calculated taking account of the ellip- 
ticity of the bulge and the inclination of the disc. Note that 
the Einstein radius does not depend strongly on whether a 
disc is included or not: there is no systematic dependence of 
the Einstein radius on the presence of a disc. 



4 LENSING PROPERTIES 

The lensing properties of discs in spiral galaxies have been 
studied by Mailer, Flores & Primack (1997), Moller & Blain 
(1998) and Blain, Moller, & Mailer (1999). In those studies 
it was found that thin exponential discs can influence the 
statistical lensing properties even if they only contain 10 per 
cent of the total mass provided that they are inclined by 
more than about 70 deg to the line of sight. In particular, the 
cross section for image geometries with 4 or more magnified 
images, and the cross section for high magnifications of 10 
or more is increased considerably. These results are generic 
and hold qualitatively for discs in early-type galaxies. How- 
ever, the scale lengths and masses of the disc components in 
early-type galaxies are likely to be smaller, relative to the 
bulge component, and so the effect of the disc on their statis- 
tical lensing properties is expected to be less. In this section, 
we describe the general and statistical lensing properties of 
discy early-type galaxies. 



4.1 General lensing properties 

The gravitational potential of a mass distribution between 
an observer at z — and a source at z — z s deflects the light 
of the source, at position /3 on the source plane, so that an 
image is observed at position 9 on the lens plane, where 



c 2 Dos 



47tGDol£>ls 



(10) 



(8) 



The deflection angle of a light ray a that passes the lens at 
an angular position 9 relative to the lens centre, is given by 



a{9) 



Dos 



(9) 



where Dos and Dls are the angular diameter distances from 
observer to source and lens to source respectively. The di- 
mensionless quantity k is defined as k — E/S c where 



is the critical surface mass density. These equations show 
explicitly why thin disc components may be important, even 
if they contain only a small fraction of the total mass: it is 
the projected surface mass density that enters equation (9) 
and this can be very high for thin edge-on discs. 

The numerical ray-tracing routines that we used to ob- 
tain the results in this paper are described in detail in Moller 
& Blain (1998, 2001). The routines were used to generate 
magnification maps on the source plane and determine the 
cross sections for high magnifications and multiple imaging. 

To illustrate the possible effect of disc components in a 
lens system, Fig. 4 shows the images of a background point 
source at a redshift of z s = 1.0 that is lensed by a fore- 
ground elliptical galaxy at redshift z\ = 0.3. The images 
are shown for the bulge-plus-disc and pure-bulge model dis- 
cussed in § 2. The source position is adjusted slightly, by 
less than 0.01 arcsec, so that both the total magnification of 
the source and the image geometries are very similar. Note 
the change in the relative fluxes of the images. The insets 
show the magnification maps on the source plane, fj,(9) for 
both cases; the disc increases the area enclosed by the caus- 
tic only slightly. Comparing the magnification maps on the 
source plane as shown in the insets of Fig. 4 with those of 
disc lenses shown in Fig. 3 in Moller & Blain (1998), makes 
clear that the effect of a disc component in early-type lensing 
galaxies is much less pronounced. However, there are differ- 
ences in the magnification ratios which we will investigate 
more in §5. 

The importance of discs can also be seen when the rel- 
ative shifts of the image positions for a source at a fixed 
position is considered. The shift gives some indication of 
how significant a given modification of a lens model is. We 
show the image geometries for small, extended sources for a 
number of our lens models in Fig. 5. There is a small, but 
measurable change in the image shapes and positions for 
E/S0 lenses of the order of 0.05 arcsec. In the case of later- 
type lenses, there is a much more pronounced change; even 
Sa lens galaxies significantly change the image geometries. 
Note that there is also a change in the overall scale set by the 
Einstein radius, of 5 — 10 per cent. This is consistent with the 
change of mass projected within the Einstein radius shown 
in Fig. 3. However, this change of scale only affects the mass 
normalization and does not reflect an observable effect of 
the disc component. 



4.2 Statistical lensing properties 

First we determine cross sections for multiple images and 
high magnifications for ellipticals that include an inclined, 
thin, exponential disc. 

In Fig. 6 we show the cross section ratios for bulge and 
bulge-plus-disc models as a function of magnification, aver- 
aged over all the galaxies of a certain morphological classifi- 
cation in the sample. Except for late-type (Sa, Sab and Sb) 
galaxies, there is no significant increase in the cross sections 
for magnifications of up to 100. Even though a few early- 
type galaxies show a clear signature of an enhancement by 
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Table 1. Elliptical galaxies in CL 1358+62. The values for < 7h >, j*h and < /hd > are taken directly from Kelson ct al. (2000). The disc 
scale length r<j is related to the half-light radius as given in Kelson et al. by r<j = 1.688rhd- in a cosmology with £7a = 0-7, !)m = 0.3 and 
#0 = SOkms" 1 Mpc -1 , 1 arcscc corresponds to 6.6 kpc at the cluster redshift of z = 0.33. A larger table that includes all the galaxies 
in CL 1358+62 investigated by Kelson et al. (2000), is given in the Appendix. 



Bulge only — > < Bulge and disc 
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Type 
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rh 






< Ih > 


rh 




< hd > 






Re 


fe 


ft 


212 


E 


20.780 


0.683 


0.488 


1.220 


21.170 


0.781 


0.340 


18.030 


0.031 


6.138 


1.160 


0.808 


0.960 


242 


E 


21.920 


1.529 


0.171 


1.360 


22.420 


1.825 


0.108 


20.720 


0.156 


0.515 


1.270 


0.806 


0.950 


256 


E 


20.880 


1.380 


0.445 


2.330 


21.200 


1.551 


0.331 


18.000 


0.053 


6.310 


2.260 


0.836 


0.967 


303 


E 


20.590 


0.638 


0.581 


1.270 


21.110 


0.774 


0.360 


18.890 


0.054 


2.780 


1.210 


0.787 


0.947 


360 


E 


20.220 


0.341 


0.817 


2.176 


21.660 


0.594 


0.217 


17.620 


0.030 


8.954 


1.000 


0.565 


0.867 


375 


E 


22.860 


4.979 


0.072 


2.650 


23.010 


5.267 


0.063 


19.220 


0.096 


2.051 


2.870 


0.800 


0.984 


409 


E 


21.270 


0.498 


0.310 


2.176 


21.490 


0.542 


0.254 


22.600 


0.132 


0.091 


0.660 


0.959 


0.969 


412 


E 


21.390 


0.767 


0.278 


2.176 


22.320 


1.091 


0.118 


17.660 


0.033 


8.630 


0.870 


0.502 


0.908 


531 


E 


21.260 


1.549 


0.313 


2.090 


21.710 


1.830 


0.207 


19.300 


0.108 


1.905 


1.980 


0.772 


0.954 


534 


E 


21.210 


0.620 


0.328 


0.860 


22.170 


0.878 


0.136 


16.680 


0.019 


1.281 


0.770 


0.489 


0.901 


536 


E 


21.180 


1.266 


0.337 


1.790 


21.520 


1.433 


0.247 


18.570 


0.057 


3.733 


1.720 


0.807 


0.965 



the disc component in the magnification map on the source 
plane, the average effect is small. 

Rusin & Tegmark (2001) showed that the fraction of 4 
image systems found in the CLASS survey is inconsistent 
with the predictions from standard lens models. They found 
that under the assumption of dark matter halos that are 
not significantly more flattened than the visible mass, about 
25 per cent of all CLASS lenses should be quads. However, 
the observed fraction is about 60 per cent. We now investi- 
gate whether disc structures in the early-type lensing galax- 
ies could explain this discrepancy. Fig. 7 shows the expected 
fraction of quad systems as a function of the total magni- 
fication threshold A for the different morphological types. 
There is no significant increase in the cross section for the 
formation of a quad lens for the subset of early-type galaxies; 
this cross section increases significantly only for the late- type 
galaxies. Many observed quad lens systems are very symmet- 
rical. This suggests that the sources do not lie very close to 
a cusp (B1422+231 is an exception), and so that the sources 
of the majority of known quad lens systems are likely to be 
magnified by only moderate amounts p ~ 10 — 20. Hence, 
discs in ellipticals cannot explain the discrepancy between 
the predicted and observed quad fraction. However, our cal- 
culations are conservative, as they do not take into account 
the proper statistical range of lens parameters, source lu- 
minosity function, the selection function of the survey and 
the magnification bias. Rusin & Tegmark (2001) show that 
such factors can increase the predicted fraction of quads at 
modest magnifications by up to a factor of ~ 2. 

Statistical lensing constraints obtained from CLASS by 
Chae (2002) suggest that the observed fraction of quad 
lenses is consistent with an average lens ellipticity e ~ 0.4, 
and thus that many SO and later lenses must be present to 
match this ellipticity. Our results in Fig. 7 suggest that if 
the majority of lens galaxies are indeed later than SO, then 
the predicted fraction of quads is consistent with the obser- 
vations. 

In summary, these results indicate that the presence of 
discs with constant mass-to-light ratio in ellipticals, is un- 
likely to affect the population statistics of strong lens sys- 
tems. However, inspecting the magnification maps on the 
source plane shows that discs might have a significant effect 



for individual lenses (Fig. 5) . To investigate this further we 
determined how much the magnification ratios of individual 
lens systems are affected when a disc component is included 
in the lens model. 



5 MAGNIFICATION RATIOS 

Attempts to model a number of known quad lens systems, 
like B1422+231 and PG1115+080 have shown that smooth 
power-law ellipsoidal models cannot fit both the observed 
image positions and magnification ratios. This has been 
generally interpreted as evidence for substructure in their 
CDM haloes (Metcalf & Madau 2001; Bradac et al. 2002; 
Chiba 2002; Dalai & Kochanek 2002; Metcalf & Zhao 2002). 
The observational detection of substructure in dark matter 
haloes would be of great importance for the standard model 
of structure formation, which predicts such CDM substruc- 
ture on galactic and sub-galactic scales (Moore et al. 1999). 
If the presence of CDM halo substructure provides the only 
possible explanation for the observed magnification ratios, 
this would provide strong observational evidence in support 
of the current CDM model. We now investigate whether 
discs in early-type lens galaxies provide an alternative ex- 
planation. 



5.1 Causes for unusual magnification ratios 

As pointed out by Mao & Schneider (1998), for point sources 
that are highly magnified and split into 4 magnified images, 
A, B, C and D, with B the brightest and D the faintest 
image (cf. Fig. 4 in the previous section; note that we differ 
here from the more conventional labelling, where images are 
labelled A,B,C and D in order of decreasing brightness), the 
individual image magnifications obey the relation 

_ MA + , 
Mac/b = — 1, 



(11) 



in the limit of utot — * oo. The loci of source positions at 
which /itot is infinite form the caustic line, and so equa- 
tion (11) holds for sources lying exactly on the caustic. What 
is not known, however, is how well equation (11) holds for 
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Figure 4. Lensing of a point source by an elliptical galaxy with 
and without a disc. The contours show the images of a point 
source at z s = 1 that lies behind a lens galaxy with the same 
properties as galaxy no. 242 in Table 1, convolved with a Gaus- 
sian beam of 0.1 arcsec FWHM. The contour levels are 0.5, 3, 10, 
20, and 30 in arbitrary flux units. The top panel shows the re- 
sulting images for a bulge-only model. The bottom panel shows 
the images for a bulge-plus-disc model, where the disc is inclined 
at i = 75 deg to the line of sight. The major axis of both bulge 
and disc is at 45 deg to the x-axis. The insets show the respective 
0.2 arcsec X 0.2 arcsec magnification maps on the source plane. 
The grey-scale is logarithmic and ranges from fj, ~ 2 (white) to 
fjL ~ 200 (black). The source position is marked by a star in the 
inset, and is slightly different in the two panels, to keep the image 
geometry very similar and the total source magnification roughly 
the same in both cases. 

sources that lie close to the caustic, or, for sources of small 
but finite size. In addition, it is not clear how the details of 
the lensing potential will affect this behaviour close to the 
caustics. The observed magnification ratios in systems like 
B1422+231 grossly violate equation (11). This discrepancy is 
the main reason why the observed magnification ratios are 
generally regarded as strong evidence for CDM-halo sub- 
structures with masses of the order of 10 6 - 10 r M , which 
can strongly affect the flux of individual images if located 
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Figure 5. Lensing of an extended source by galaxies with and 
without discs. The grey-scale shows the images of a small ex- 
tended source of size ~ 0.1 arcsec at z B = 1 that lies behind lens 
galaxies with the properties of the labelled galaxies in Table Al. 
The models in the right panels include a disc at an inclination 
angle of i = 75 deg: the bulge and disc components are oriented 
at 45 deg with respect to the x-axis in all cases. The left panels 
show the corresponding images for bulge-only models. 



close to the images. However, there are several other expla- 
nations: 

(i) microlensing. In many strong lens systems, microlens- 
ing variability has been detected and it is difficult to exclude 
the possibility of microlensing in the observed systems with 
certainty. 

(ii) differential extinction. There is a possibility that op- 
tical flux ratios of the images are affected by different dust 
columns along the lines of sight through the lensing galaxy. 

(iii) the source is extended. Extended radio sources may 
lead to very different magnification ratios in the observed im- 
ages, depending on the size and relative brightness of the dif- 
ferent source components. Recent simulations by Moustakas 
& Metcalf (2003) suggest that even for extended sources 
there is no change in the continuum magnification ratios 
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Figure 6. The ratio of the magnification cross sections. The 
curves show the ratio of the cross sections, o"b+d(M > A)/<r^,(fj, > 
A), for total magnifications fi > A of a point source at z B = 1 as 
a function of the threshold magnification, A. The different curves 
are calculated by averaging over all the galaxies in Table Al that 
have the same morphological type. When a disc is included, the 
inclination angle i = 75deg in all cases. 
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Figure 7. The fraction of four-image lens systems as a function of 
minimum magnification for different lens models. For each mor- 
phological class of galaxies listed in Table Al the average cross 
section ratio <r N _4/a- tot is shown, where on = 4 * s the cross section 
for the formation of 4 images and ctot is the total cross section 
for multiple imaging. The left panel shows the results for earlier 
types, the right panel for later types. In all bulgc-plus-disc mod- 
els, a disc with the parameters as listed in Table Al is included, 
with an inclination angle i = 75deg. 



for smooth lens models; only the relative contributions from 
broad and narrow lines is affected. 

(iv) off-caustic sources. If the source does not lie exactly 
on the caustic but is still magnified strongly, then the mag- 
nification ratios will not necessarily obey equation (11). 

We now investigate the validity of the last explanation in 
the context of discs in early-type lensing galaxies. 



5.2 Known lens systems with unusual 
magnification ratios 

There have been two recent studies of lens systems with 
unusual magnification ratios. In the comprehensive study 
by Metcalf & Zhao (2002) 4 systems are listed where 
the magnification ratios deviate by ~ 0.1 mag or more 
from the expected values assuming a one component 
smooth ellipsoidal lens with external shear; B1422+231, 
PG1115+080, H1413+117 and Q2237+030. The first three 
systems are quasars lensed by early-type galaxies while the 
lens in Q2237+030 is a barred spiral. Of the four systems, 
HT413+117 lies close to a cluster of galaxies which may cru- 
cially affect the lensing potential (Kneib et al. 1998; Chae & 
Turnshek 1999). Q2237+030 is lensed by a galaxy at a very 
low redshift of z — 0.04 and microlensing affects the flux ra- 
tios by up to a factor of two (Irwin et al. 1989; Wozniak et al. 
2000). PG1115+080 and B1422+231 are affected only very 
moderately if at all by microlensing (Vanderriest et al. 1986; 
Yee & Bechtold 1996), but both reside very close to, or even 
within, a galaxy group (Tonry 1998; Kundic et al. 1997a,b). 
As shown by Moller et al. (2002), the presence of the group 
can significantly alter the lensing behaviour and it is doubt- 
ful whether a simple power-law-plus-shear model, as used 
for example by Keeton, Kochanek & Seljak (1997) reflects 
the true lensing potential adequately. For B1422+231 the 
predicted and observed magnification ratios differ by up to 
0.3 mag, as compared to only 0.1 mag for PG1115+080, and 
so it is less likely that the discrepancy in B1422+231 is a re- 
sult of the details of the group mass distribution surrounding 
the lens. This points out B1422+231 as the strongest candi- 
date for dark substructures: it is studied in great detail by 
Bradac et al. (2002). 

A recent theoretical study (Keeton, Gaudi & Petters 
2002) investigated the magnification ratios of four caustic 
lens systems in more detail, using an analysis similar to that 
presented here but with a simple one-component lens model: 
B2045+265, B0712+472, RX J0911+0551 and B1422+231 
(Keeton, Gaudi & Petters 2002). RX J091 1+0551 has only 
optical and near-infrared flux measurements, and so the 
extinction corrections for this source are highly uncertain. 
B0712+472 is possibly affected by a foreground group (Fass- 
nacht & Lubin 2002). Keck spectra of B2045+265 indicate 
that it is likely to be a late- type lens (Fassnacht et al. 1999), 
and so it is unlikely that a simple one-component model is 
viable: a disc is almost certainly present. Note that only 
B1422+231 was included in the study of Metcalf & Zhao 
(2002). 

The number of lens systems that are claimed to have 
magnification ratios that can only be explained using com- 
pact CDM substructure is rather large. In some cases, mi- 
crolensing variability or differential extinction could be a 
viable alternative explanation, but not for all. 

We summarise the properties of all these lens systems 
in Table 2. Together with the image fluxes, we list the mag- 
nification ratios /i A c/B and /iabc/d = (ma + Mb + Hc)/nn- 



5.3 Magnification ratios for bulge-plus-disc 
models of early-type lens galaxies 

To answer the question whether discs in early-type lenses 
could explain some or all of the observed magnification ra- 
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Table 2. Summary of observed quad systems with unusual magnification ratios. 1 Christian, Crabtree & Waddcll 1987. The fluxes are 
normalized to image D. 2 Jackson ct al. 1998. The listed values are for the Merlin 1996 5 GHz data. 3 Fassnacht et al. 1999. The listed 
values correspond to the VLA 8.5GHz data. 4 Patnaik et al. 1999. The listed values correspond to the VLA 8.5GHz data. 5 Chae & 
Turnshek 1999. All images have been relabelled to correspond to the convention adopted in other systems. 6 Falco ct al. 1996. All images 
have been relabelled to correspond to the convention adopted in other systems. 7 Burud et al. 1998. The fluxes are normalized to image 
A. 



System 


Fa 


Fb 


F C 


Fu 


Mac/b 


^ABC/D 


PG1115+080 1 


2.49±0.03 


3.22±0.03 


0.64±0.03 


1.0 


0.98±0.1 


6.35±0.1 


B0712+472 2 


10.5±0.1m.Jy 


8.5±0.1mjy 


4.7±0.1mjy 


0.9±0.1mJy 


1.78±0.03 


26.3±2.1 


B2045+265 3 


18.4±0.3mJy 


9.42±0.3mJy 


14.8±0.3mJy 


2.41±0.3 mjy 


3.5±0.1 


17.7±2.1 


B1422+231 4 


152±2 


164±2 


81±1 


5±0.5 


1.42±0.03 


79.4±7.9 


H1413+117 5 


1.65±0.16 


1.62±0.16 


1.0 


0.67±0.1 


1.63±0.16 


6.1±0.6 


Q2237+030 6 


78±14 ^tjy 


60±14 ^tjy 


85±14 /ijy 


43±14 /ijy 


2.7±0.8 


5.2±2.4 


RX J0911+0551 7 


1.0 


0.965±0.013 


0.544±0.025 


0.458±0.004 


1.60±0.08 


5.48±0.28 



tios, wc computed all the image multiplicities, position and 
magnifications for all sources on a 400 x 400 pixel grid for 
our particular representative model galaxy (galaxy no. 242). 
We performed the calculations both on a low-resolution grid 
that contained the whole caustic structure and on a smaller, 
higher resolution grid, that contained one of the cusps along 
the major axis of the lens galaxy. For both grids, we se- 
lected all sources with quad images, yielding a total sample 
of about 40000 simulated image systems. 

Fig. 8 shows a scatter plot of the ratio of the sum of 
the magnifications of the outer two images to the central 
image /iac/b, versus the ratio of the total magnification of 
the brightest three images, A-C to the magnification of the 
faintest image D, /iabc/d- The points in the upper panel are 
for the pure-bulge model, whereas the points in the lower 
panel are for the bulge-plus-disc model. For large magnifica- 
tion ratios there is a significant difference. All the images of 
sources with high total magnifications, corresponding to a 
high /iabc/d ratio, scatter around ^ac/b ~ 1 for the bulge- 
only models. For large the scatter of a few per cent equals 
that expected from numerical uncertainties in the simula- 
tions, due to small deviations in the deflection angle that 
lead to large deviations in the source position and images 
occasionally being missed or associated with a wrong source 
position. Note that there are two distinct and well defined 
'populations' of images, corresponding to images in the two 
different types of cusp regions (see Fig. 9). The two panels 
of Fig. 9 show quantitatively how much the image magni- 
fication ratios deviate from equation (11) as a function of 
source position, both for the bulge-plus-disc model (upper 
panel) and the bulge-only model (lower panel). These maps 
are similar overall, but differ significantly in the two cusps 
that are aligned with the disc: when a disc is present, im- 
ages of sources that lie very close to these cusps have mag- 
nification ratios that are inconsistent with the predictions 
of equation (11). The images that are produced by sources 
that lie close to these two cusps form the second population 
of images in the bottom panel of Fig. 8. 

As an example, Fig. 4 illustrates the behaviour for 
sources in the lower-left cusp of our model. By comparing 
the relative fluxes of the images in the top and bottom pan- 
els of Fig. 4, one can see by eye that the sum of the fluxes 
of the images A and C is roughly equal to the flux of image 
B in the bulge only model, but not in the bulge-plus-disc 
model. The relative positions and flux ratios of the images 
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Figure 8. Scatter diagrams of magnification ratios for quad im- 
ages of point sources lensed by a model galaxy with parameters 
matching galaxy no. 242 in Table 1. The markers in the upper 
panel show the magnification ratios between images produced by 
the bulge-only model, whereas the markers in the bottom panel 
show the results for the bulge-plus-disc model. For total magnifi- 
cations below 20 only 10 per cent of all simulated image systems 
are shown. Note that the images in the bulgc-plus-disc case fall 
into two regions, whereas all the images produced by pure-bulge 
models obey equation (11) in the limit of MABC/D — ► 1- The large 
markers indicate the magnification ratios for the observed systems 
with the la error, as listed in Table 2. A small scatter of about 
10 per cent is introduced at high fi due to numerical uncertainties. 



A,B and C in the lower panel of that figure are very close 
to those observed in B1422+231. 

Figs 8 and 9 show the difference between the bulge and 
bulge-plus-disc models. In the bulge-plus-disc model a set of 
sources exists with an image magnification ratio /iac /b very 
different from one, even for high magnifications of /iabc/d ~ 
50. All these sources lie inside one of the two cusps along the 
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Figure 9. Maps of the magnification ratio [(/^a + McVmb — 1] — 
as a function of source position. The upper panel is for the bulge- 
only model, the lower panel for the bulge-plus-disc model. Large 
values (dark shades) indicate that the magnification ratios of the 
images of a source at that position obey equation (11). Small val- 
ues (light shades) indicate that the image magnification ratios 
deviate from that relation. Note that Fig. 4 shows an example 
configuration for sources lying in the bottom left cusp. 



major axis. Note that for all those images, the faintest image 
D lies at a large distance of ~ 1.2 arcsec from the centre of 
the lensing galaxy and is not demagnified. This reflects the 
high total magnification of a source which lies close to the 
caustic line. 

Some unusual observed magnification ratios of some ob- 
served lens systems are indicated by the large symbols in 
Fig. 8. Error bars indicate the 1-a flux measurement errors 
(Table 2). Only H1413+117 and RXJ0911+0551 lie in re- 
gions of the diagram that are consistent with a bulge-only 
model. 

However, the bulge-plus-disc model is clearly consistent 
with the magnification ratios for PG1115+080, H1413+117, 
RX J0911+0551, B0712+472, Q2237+030 and B1422+231. 
B2045+265 is the only system for which our model cannot 
reproduce the observed magnification ratios. This source is 
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Figure 10. Scatter diagrams of magnification ratios for quad 
images of point sources lensed by a model galaxy with parameters 
as in the Fig. 9, but with a heavier disc that has twice the size. 
The symbols are in the previous plot. 



lensed by a spiral galaxy, and so our early-type lens model 
with a very modest disc would not be expected to provide a 
good fit. 

In Fig. 10 we show magnification ratios for about 40000 
simulated image systems (without using a second higher- 
resolution grid) for a bulge-plus-disc model with a disc that 
has a scale-length twice that of the previous bulge-plus-disc 
model. All the other parameters are unchanged and corre- 
spond to galaxy no. 242. The effect of increasing the disc 
mass is to substantially increase the hac/b ratio for image 
systems with magnifications Mabc/d ~ 20. The magnifica- 
tion ratio of B2045+265 is now consistent with the predicted 
value. 

We conclude that the presence of an inclined disc com- 
ponent containing of order 5 per cent of the total mass can 
describe all of the known early-type lens systems with un- 
usual magnification ratios without requiring compact sub- 
structure in the lens halo. 



5.4 Effect of disc size, disc inclination, bulge 
ellipticity and a dark halo 

We have used a particular model galaxy (§ 2) as a lens. In 
order to make sure that our results apply in general to bulge- 
plus-disc models of early-type galaxies and are not specific to 
our particular choice of parameters we investigated a range 
of other models. The full parameter space for bulge-plus-disc 
models is large and a detailed study is beyond the scope of 
this paper. We investigated the dependence on several pa- 
rameters qualitatively. We find that, qualitatively, any disc 
that on its own produces a caustic of similar or larger size 
than that produced by the bulge is effective at modifying 
the image properties. In practice that means any disc that 
is inclined by more than 60 deg with a scale-length of the 
order of the size of the caustic of the corresponding bulge- 
only model. For nearly spherical bulges, discs with smaller 
scale- lengths are also effective, and naturally produce mag- 
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nification ratios similar to those observed for B1422+231. 
For disc sizes below 0.2 arcsec, the disc inclination must be 
increased beyond i = 80 deg to obtain magnification ratios 
consistent with Bf422+23f. Even an edge-on disc smaller 
than 0.1 arcsec cannot reproduce the magnification ratios of 
B1422+231. There is also the possibility of misalignment 
between the positions of the disc and the bulge/halo. Such 
misalignments are to be expected as a result of recent merger 
events and may lead to complicated and unusual lensing con- 
figurations and provide more extreme changes in the magni- 
fication ratios, as discussed by Quadri, Moller & Natarajan 
(2002). 

We have also investigated models with a CDM halo. 
Adding a singular CDM halo, containing 2/3 of the total 
mass within an Einstein radius, and adjusting the mass nor- 
malization of the bulge and disc so as to keep the Einstein 
radius constant, does not significantly alter the effect of the 
disc. It is purely the presence of a thin disc structure in 
combination with a more spherical mass component that 
produces the change in magnification ratios. 

5.5 Expected number of affected systems 

Moller & Blain (1998) have shown that discs in spiral galax- 
ies are important for the statistical effects on lensing cross 
sections only at inclination angles of > 60 deg. Assuming 
that a fraction of 50 — 100 per cent of all early-type lens 
galaxies do contain discs, 15 — 35 per cent will contain discs 
with inclination angles of > 60 deg. A conservative estimate 
of the fraction of early-type lens systems in which discs may 
be important is thus ~ 10 — 30 per cent. This fraction is 
consistent with the number of CLASS quad lens systems 
with deviant magnification ratios: of the 7 CLASS quad 
lenses, at least two early-type lens system (B1422+231 and 
B0712+472) exhibit magnification ratios inconsistent with 
simple power-law plus external-shear models. 

An additional factor to consider is magnification bias. 
Depending on the bulge-to-disc ratio, the disc may lead to 
higher cross sections for strong magnifications, as shown in 
Fig. 6. This could bias a flux-limited sample towards lens 
galaxies that contain discs. However, as shown in Fig. 6, the 
cross section is not increased significantly for E/SO/Sa galax- 
ies, and so the magnification bias is likely to be small. 

5.6 Distinguishing the effect of discs and CDM 
substructure 

The results described in this section demonstrate that discs 
in early-type lensing galaxies may explain some of the ob- 
served magnification ratios that have hitherto only been ex- 
plained assuming CDM substructure. There are some key 
observational differences between these two scenarios. 

First, discs in early-type lens galaxies are, in principle, 
detectable by bulge-plus-disc light profile decompositions as 
discussed in § 2. However, the light of the background quasar 
makes accurate photometry of lens galaxies extremely dif- 
ficult in practice. Treu & Koopmans (2002) obtained kine- 
matic information of a lensing galaxy from spectroscopic ob- 
servations, and demonstrated that it is thus possible to ob- 
tain additional constraints on its mass distribution. It might 
be possible to constrain possible disc structures in lensing 
galaxies in this way. 



Secondly, disc structures significantly affect only the 
magnification ratios of sources that lie within one of the two 
cusps along the discs. Assuming that discs are oriented along 
the major axis of the bulge ellipticity, this means that disc 
structures provide an explanation of discrepant magnifica- 
tion ratios only for quad image systems with three strongly 
magnified images along a line perpendicular to the major 
axis of the lens galaxy, whereas CDM substructure could 
modify any image configuration. Note that this is not true 
for late- type lenses with extended inclined discs, in which 
the caustic structures are completely dominated by the disc, 
and the disc strongly affects the magnification ratios of quad 
systems produced by sources in any one of the four cusps. 

We stress that these results do not prove the existence 
of discs in early-type lensing galaxies, nor do we propose 
that such disc structures provide the only explanation for 
all the observed magnification ratios. For all systems, the 
observed magnification ratios could be caused by CDM sub- 
structures. We have not carried out detailed modelling of any 
of the observed systems, and our analysis does not include 
joint fitting to both magnification ratios and image posi- 
tions. However, in contrast to Keeton et al. (2002) we find 
that magnification ratios alone do not provide conclusive ev- 
idence for the presence of CDM halo substructure: realistic 
disc structures can explain the observed magnification ra- 
tios equally well. Detailed case-by-case lens modelling using 
bulge-plus-disc models is now needed to verify that early- 
type galaxies with disc structures are indeed viable models 
for all or at least some of the observed lens systems. The fact 
that discs affect image positions very slightly, even for a fixed 
source position, indicate that detailed bulge-plus-disc mod- 
els are likely to be successful in fitting magnification ratios 
and image positions. 



6 EFFECT ON HUBBLE PARAMETER 
ESTIMATION 

Strong gravitational lensing of variable background quasars 
provides a unique method to determine the Hubble constant. 
Brightness fluctuations are observed at different times at 
each image as the path lengths and gravitational time delays 
to each image depend on its relative position and the local 
lensing potential. If the image positions and lensing poten- 
tial are known, then Ho can be determined, as At oc 1/Ho- 
Gravitational and geometric time delays are often of roughly 
equal magnitude, and so an accurate knowledge of the lens- 
ing potential is crucial to deduce an accurate value of Ho- 

Lens systems with measured time delays are nearly 
always associated with early-type lens galaxies, although 
B1600+434 (Koopmans et al. 2000) and B0218+357 (Wuck- 
nitz, priv. comm.) are notable exceptions. These systems 
have been modelled using bulge-only models, with a con- 
stant external shear in some cases. 

We now investigate how a disc component is likely to 
affect time delays for the bulge and bulge-plus-disc model 
of galaxy no. 242. The time delay between two images at 9\ 
and 02 of a source at (3 is 



where Dql is the angular diameter distances between ob- 
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Figure 11. Time delay as function of source position. The two 
maps show the time delays between images B and D (cf. Fig. 4) 
for all quad geometries for model galaxy no. 242. The left panel 
shows the results for the bulge-only model and the right panel the 
results for the bulge-plus-disc model. The boxes in the left panel 
mark the different cusp regions, as discussed in the text. 

server and lens. The distances and thus At are proportional 
to Hq 1 . The lens redshift is Z\ and the potentials due to 
bulge and disc components are iph and Vd respectively. In 
§ 4 we pointed out that the shift in image positions between 
bulge-plus-disc and bulge-only models for E and E/SO type 
lensing galaxies is expected to be very small. 

We are concerned here with the effects of modelling an 
early-type lens that contains a disc with a bulge-only model. 
Since the lens models are obtained from fits to the image 
positions, the image positions remain the same irrespective 
of the model used. Therefore, the geometrical time delay of 
a given pair of images is not affected by a disc component; 
any change in the time delay is due to the change in the 
potential at the image positions. 

Note that compact CDM substructure does not affect 
time-delays measurably. 



6.1 Four-image/quad systems 

About half of the known lens systems with measured time 
delays are quad systems. Since the quad caustic structure is 
the region in the source plane that is most strongly affected 
by the presence of a disc, the greatest effect of a disc on 
time delays is expected in these systems. In most cases a 
single time delay is measured, usually between the image 
pair with the largest separation. Using our ray-tracing code 
in combination with equation (12) we calculated the time 
delays between images B and D for all quad systems of the 
model galaxy described in § 2. Fig. 11 shows the time delays 
as a function of source position for bulge-plus-disc and bulge- 
only models. The difference is on the order of a few tens of 
per cent, and is largest in the cusps lying along the disc (the 
regions marked I). Especially in the cusp regions close to 
the caustic, where the image magnification is greatest, the 
time delays for the bulge-plus-disc model are 20 — 30 per cent 
larger than for the bulge-only model. 

The time-delay statistics for regions I and II are shown 
in Fig. 12. For each region in Fig. 11 the time delays are 
divided into 50 bins and the histograms are normalized 
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Figure 12. Histogram of time delays. This plot shows the ex- 
pected relative fraction of systems as a function of timedelay be- 
tween images B and D. The left panel shows the statistics for 
the regions marked I in Fig. 11 and the right panel the statis- 
tics for the regions marked II. The histograms are shown for our 
bulge-plus-disc and bulge-only lens galaxy model (see §2). Note 
that the time delay is smaller for bulge-only systems, due to the 
shallower potential. 
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Figure 13. Histogram of the separations of images. This plot 
shows the expected relative fraction of systems as a function of 
separation. Note that the separation is larger for bulge-only sys- 
tems, due to the higher mass in the central region. 



to unity. Comparing the histograms in the right panel of 
Fig. 12, there is no significant difference between the bulge- 
plus-disc and bulge-only model time-delay statistics in the 
cusps away from the disc (the regions marked II on Fig. 11). 
However, the histograms in the left panel of Fig. 12 show 
that the time delays in the cusps lying along the disc are in- 
creased by ~ 20 per cent for the bulge-plus-disc model with 
respect to the bulge-only model. The overall distribution for 
the bulge only model is shifted towards shorter time delays 
with respect to the bulge-plus-disc models. Assuming that 
an elliptical lens galaxy that contains a disc is modelled us- 
ing a bulge-only model, this suggests that the Hubble con- 
stant will be underestimated systematically. 

In Fig. 13 we show explicitly that disc components al- 
ter the individual image positions only very slightly for our 
model galaxy no. 242. The figure shows the distribution of 
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Figure 14. Histogram of time delays for 2-image systems. This 
plot shows the expected relative fraction of systems as a function 
of time-delay between the two images in 2-image lens systems. 
The histograms are shown for our bulge-plus-disc and bulge-only 
lens galaxy model as described in § 2. 



image separations for both models. The images are system- 
atically further apart in the bulge-only models, suggesting a 
larger mass inside the Einstein radius and thus a larger pre- 
dicted time delay, if the radial profile of the potential is unaf- 
fected. In observed lens systems, the lens would be modelled 
to fit the image configurations, and so a bulge-only model 
fitted to the observed image separations would contain less 
mass within the Einstein radius than the bulge-only model 
with the time delay distribution shown in Fig. 12. This would 
shift the time-delay distributions shown in Fig. 12 for the 
bulge-only model to even lower time delays. The only pos- 
sibility both to fit the image separations and obtain similar 
distribution of time delays is by changing the gradient of 
the potential between the images. For point sources, the 
only other lensing property that is sensitive to the shape 
of the potential is the magnification ratio between the im- 
ages. Thus, if magnification ratios are not used to constrain 
lens models, then bulge-only models might fit the image po- 
sitions of a lens that actually consists of a bulge-plus-disc 
lens, yielding an incorrect mass model and an Ho value that 
is systematically too low by up to 25 per cent. 



6.2 Two-image systems 

Quad lens systems are generally regarded as being more 
useful than 2-image systems for Ho determinations because 
more constraints on the lensing potential are available. How- 
ever, many 2-image lens systems have been used for time de- 
lay measurements. In these cases, the source lies outside the 
diamond-shaped caustic, where magnifications /u ~ 2 — 10. 
These images are expected to be affected only weakly by a 
disc component. 

We have calculated the time delays for 2-image systems 
and show the corresponding statistics in Fig. 14. Comparison 
with Fig. 12 shows that the effect of a disc component on the 
time delay is much smaller for 2-image systems than for quad 
systems. A number of 2-image systems with measured time 
delays, like B0218+357 (Biggs et al. 1999) have possible disc 
structures, but even if these are present, our results indicate 



that the the effect on estimates of Ho is small. Note that 
recent observations (Wucknitz, private comm.) indicate that 
the lens in B0218+357 is a spiral galaxy. Massive discs, as 
expected in spiral lenses, may affect the time delays much 
more strongly. 

6.3 Parameter dependence 

The results in this and the preceding section show that any 
inclined discs within elliptical lens galaxies affect sources 
in the cusps lying along the disc most strongly, whereas the 
effect on the cusps lying away from the disc is small. We have 
not checked the dependence of the time-delay statistics on 
the disc parameters in detail, but the dependence is similar 
to that discussed in § 5.4 since magnification ratios and time 
delays are both sensitive on changes in the cusps. The results 
we obtained in that section also hold for the time delays. 
Disc structures are important if they are inclined at i > 
60deg to the line of sight and are massive/extended enough 
to produce a caustic of a size that is of the order of the 
caustic produced by the bulge. That is typically the case 
for discs larger than rhd > 0.6kpc and with masses of at 
least ~ 5 per cent of the bulge, for which time delays will 
be increased by a few tens of per cent relative to bulge-only 
models. 



6.4 The effect on estimates on Ho 

Our approach here only gives a rough indication of how de- 
terminations of Ho are affected by disc structures. To be 
any more quantitative about the uncertainties that arise it 
is necessary to perform detailed modelling of known systems 
using bulge and bulge-plus-disc models, which is difficult and 
time-consuming. However, our results offer a possible reso- 
lution as to why there is a significant discrepancy between 
determinations of Ho from CMB data (Lewis & Bridle 2002), 
the HST Key-project (Freedman et al. 2001) and gravita- 
tional lensing (Kochanek 2003). 

The effect of disc structures on time-delays is much 
more significant for quad systems than for 2-image sys- 
tems. In addition, we expect that there will be a correla- 
tion between low- .Ho results and unusual magnification ra- 
tios if discs are present. This seems to be the case: Biggs 
ct al. (1999) derive H = 69±H kms -1 Mpc -1 for the 2- 
image system B0218+357, whereas Impey et al. (1998) ob- 
tain H = 44±4kms~ 1 Mpc -1 for PG1115+080, a quad sys- 
tem with unusual magnification ratios. The errors on the Ho 
estimate by Biggs et al. (1999) have recently been disputed 
by Lehar et al.(2000) and Wucknitz, Biggs & Browne (MN- 
RAS, submitted). More recent single-component models of 
PG1115+080 by Treu & Koopmans (2002) give a higher 
value of Ho = 59l^ 2 ±3 kms -1 Mpc -1 . The larger errors are 
statistical, and the smaller error of ±3 is an estimate of the 
systematic error. No other quad systems with measured time 
delays have unusual magnification ratios, and the values of 
Ho determined from these systems are generally higher. 



7 CONCLUSIONS 

Gravitational lensing plays a role of ever increasing impor- 
tance in determining cosmological parameters, notably Ho, 
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and constraining the dark matter distribution on a variety of 
scales. More than 50 strong gravitational lens systems are 
known. Generally, it is not possible to constrain the mass 
distribution of the lensing galaxy uniquely, due to the de- 
generacies present. Whether this uncertainty in models has 
serious implications for cosmological parameter estimation 
from strong lensing is still unclear. 

The majority of all lensing galaxies have been classified 
as massive early-type galaxies. Reliable photometric deter- 
mination of the presence of a disc component is difficult due 
to the usually very bright lensed quasar images. Therefore, 
discs similar to those discussed in this paper would not gen- 
erally have been identified in existing images. However, ob- 
servations of early-type galaxies in the local universe and in 
the moderate-redshift cluster CL 1358+62 suggest strongly 
that a large fraction of early-type galaxies contain disc-like 
structures. 

We have investigated how the lensing properties of 
early-type galaxies would change if a thin disc component 
that contains about 5 per cent of the mass is included: 

(i) If early-type galaxies contain disc components, the sta- 
tistical lensing properties of early-type lensing galaxies are 
affected by less than about 5 per cent. Only if the major- 
ity of lensing galaxies are of SO/Sa or later types would the 
expected relative number of quad to 2-image lens systems 
increase by a factor of 2 or more. 

(ii) A disc component in early-type lensing galaxies af- 
fects significantly the expected magnification ratios of highly 
magnified quad systems. If a disc component is present in 
the lens and inclined at more than ~ 70 deg to the line of 
sight, the resulting image configurations and magnification 
ratios can resemble that of systems like B1422+231. 

(iii) Time delays may be affected significantly by the pres- 
ence of a disc component in individual early-type lens galax- 
ies. Bulge-only lens models used to fit bulge-plus-disc lenses 
fit lensed image positions well, but would yield a value for 
Hq that is systematically low by about 25 per cent. 

The presence of inclined exponential discs in early-type 
lensing galaxies affect the lensing properties of quad lens sys- 
tems significantly. Assuming that the majority of all early- 
type lens galaxies contain discs, then both the time delays 
and magnification ratios are expected to be affected by disc 
components in 10 — 30 per cent of all quad lens systems. One 
or two of the quad lenses identified by CLASS is expected to 
be affected (consistent with the properties of B1422+231). 
Bulge-plus-disc models provide an alternative explanation 
to the observed magnification ratios and lead to systemat- 
ically longer time-delays. In order to obtain more stringent 
constraints on the presence of CDM substructure in the lens 
and on the effect of discs on Ho measurements, detailed lens 
modelling of known lens systems will be necessary. Fitting 
models that include exponential discs to real lens systems is 
much more computationally intensive than the calculations 
performed here, but may be possible. 

We have discussed the first order effect of discs in early- 
type galaxies on observed lensing properties. Once a large 
number of lens systems is known, it might also be possible to 
use lensing to constrain the masses and parameters of those 
systems. Determining the abundance and properties of discs 
in ellipticals in more detail, and at different redshifts, using 



gravitational lensing would shed new light on the details of 
the formation of discs in early-type galaxies. 
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APPENDIX A: GALAXIES IN CL 1358+62 

All the galaxies for which Kelson et al. (2000) provides sur- 
face brightness profile information are employed in our cal- 
culation of certain statistical properties, like the high magni- 
fication cross sections in §4.2. For completeness, we present 
here the surface brightness profile parameters for all the 
galaxies. 
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Table Al. Galaxies in CL 1358+62. The values for < /j, >, and < I^d > are taken directly from Kelson et al. (2000). The disc scale 
length r<j is related to the half-light radius as given in Kelson et al. by = 1.688fhd- 111 a cosmology with Q\ = 0.7, Qm = 0.3 and 
H = SOkms" 1 Mpc -1 , 1 arcsec corresponds to 6.6 kpc at a redshift of z = 0.33. 



